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ABSTRACT 

We show that galaxy clusters conform to a tight relation between X-ray isophotal size Rj and emission 
weighted intracluster medium (ICM) temperature {Tx) ■ The best fit relation for 41 members of an X- 
ray flux limited cluster sample is: logRj = (0.93 ± OTl) log {{Tx) /6keV) - (0.08 ± 0.01); intrinsic 
scatter in size about the relation is 15%, and for 30 clusters with (Tx) > 4 keV the scatter is reduced to 
10%. The existence of the size-temperature (ST) relation indicates the ICM structure is a well behaved 
function of (Tx). We use an ensemble of gasdynamic simulations to demonstrate that a cluster population 
experiencing present epoch growth nonetheless conforms to an ST relation with scatter similar to that 
observed; the simulations also exhibit a tight relation between Myir and (Tx) , providing the suggestion 
that a similar relation holds for observed clusters. We use the scatter in Rj to estimate limits on the 
rms variation in ICM mass fraction SficM at constant (Tx) : 5f / ficM < 22% (< 14% for clusters 
with (Tx) > 4 keV). It appears that a mechanism like feedback from galaxy winds, which introduces 
systematic structural changes in the ICM, is required to reproduce the observed slope of the ST relation. 

Subject headings: galaxies: clusters: general — intergalactic medium — X-rays: galaxies 



1. INTRODUCTION 



Ongoing cluster growth at the present epoch i s an ob 
servational fact (e.g 



Fabricant et al. 1986, 



off fc Zaritoky 1005 | , pcnry fc Briel 1996 | , pMohr, CcUer & 



1989, Zablud 



Wogncr 199C )4- largo sFudics employing varied tochniqu cE 



indicate that 30-70% of clusters e x hibit evidence for re- 
cent mergers (peller fc Beers 1982|, ^ressler fc Shectmar^ 



19881 pones fc Forman 1992| , |Mohr et al. 199^ ). A theoret- 
ical framework for understanding this growth has emerged, 
and structure formation simulations within a range of fa- 
vored cosmological models produce clusters which exhibit 
merger si gnatures similar in detail and frequen cy to those 
observed ( |Mohr et al. 1995| , puote fc Xu 19971 ). 

Yet, recent numerical work indicates that even though 
clusters are growing at the present epoch, they are, on 
average, regular objects with temperatures and masses 



tightly correlated through the virial theorem (Evrard,| 
Metzler fc Navarro 1996 , Schindler 1996| , Roettiger, Burn;^ 



fc Loken 1996). Direct observational tests of these results 
are difficult, but analyzing relations between ICM struc- 
tural parameters and temperature provides an alternative 
means of testing cluster regularity. The well known X-ray 
luminosity-temperature (Lx — (Tx) ) relation (e.g. smith 
Mushotzky fc Scrlemitsos 1979|,|MitclicU et al. 1979|, |Davic 



et al. 1993| , jMushotzky fc Scharf 199"^ is one observationa 
indicator of cluster regularity, but the scatter around the 
relation is quite large. The scatter is largely due to vary- 
ing degrees of excess core emission associated with cooling 
ffows (Fabian et al 1994). 

Here we report a tight relation between cluster isopho- 
tal size Ri and mean, emission weighted temperature (Tx) 
that we uncovered in an ongoing morphological analysis of 
a large. X-ray flux limited cluster sample. This ST relation 
for observed clusters indicates a high degree of regularity, 
suggests the existence of a mass-temperature correlation, 
and provides new constraints on theories of the interac- 



tions between galaxies and the ICM. This newly discov- 
ered ST relation may ultimately lead to more accurate es- 
timates of cluster baryon fractions and tight er constraints 



on the cosmic den sity parameter VLq (e.g. White et al 
1993| , [Evrard 1997| ) 



Section 2 contains a description of the data reduction, 
analysis and the observed ST relation. We then use an 
ensemble of numerical simulations to examine the ST re- 
lation in a cluster population experiencing present epoch 
growth (§3). In Section 4, we discuss sources of scatter 
about the relation, estimate upper limits on the variation 
in ICM mass fraction ficM , and discuss the slope of the 
ST relation. Section 5 contains a discussion of our conclu- 



2. OBSERVATIONS 

The cluster sample contains members of the X-ray flux 
limited group of 55 clusters defined by Edge et al. (1990). 
Archival ROSAT Position Sensitive Proportional Counter 
(PSPC) images of 47 of these clusters are available online 
through the High Energy Astrophysics Science Archive Re- 
search Center (HEASARC). We reduce thes e images us- 
ing PRO S and Snowden analysis software (Snowden et 
al. 1994). The final image for each observation is the 



sum of individually flatfielded subimages from the Snow- 
den bands R4 through R7, corresponding to photon ener- 
gies 0.44-2.04 keV. We exclude time intervals with master 
veto rates higher than 220 cts/s, and exclude other high 
background time intervals (typically <5 min excluded) 
whose inclusion would degrade the detectio n significance 
of a source 10% as br ight as the background ( Pildis, Breg-| 
man fc Evrard 1995). We combine multiple images of a 



cluster using the positions of bright X-ray point sources 
where possible, or alternatively, the image header point- 
ing positions. Finally, we remove obvious point sources 
from the unsmoothed images (we do not remove possi- 
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ble point source components coincident with the extended 
emission typical of cluster cooling flows, as in NGC 1275, 
for example) and then Gaussian smooth (cr = 2 pixels). 
The pixel scale is 14". 947, and the effective angular reso- 
lution is FWHM ^ 1.5'. We correct the cluster images 
for (l + z)^ cosmological dimming, but do not make galac- 
tic absorption or K corrections (discussed in more detail 
below). 

We define the cluster size Rj using the area A/ of the 
largest region enclosed by the isophote /. 



(2-1) 



This approach deals consistently with clusters that are not 
azimuthally symmetric because of recent mergers, and pro- 
duces a tighter relation than would azimuthal averaging; 
in addition, / can be chosen so that Ri is unaffected by 
the central surface brightness excesses typical of cooling 
flows. The isophote / is the cluster surface brightness, 
and we determine the X-ray background Ib using annuli 
well outside the X~ray bright cluster region (typically 38' 
from the cluster center with ^ 8' extent). We estimate 1b 
as the (sigma clipped) mean value of the pixels within the 
annulus, and we use the width of the distribution around 
Ib to estimate the uncertainty. For the isophotes / we use, 
the uncertainty in Ri is primarily due to the background 
uncertainty. We conservatively estimate the uncertainties 
in Ri to be an ^ {Ri_^ - Rj^) /2 where /+/_ = / ± /s/lO. 
The internal estimate of the uncertainty in Ib is signifi- 
cantly smaller than 10% in all cases. 

We use published, emission weighted, mean tempera- 
tures (Tx) for each cluster. There are 42 clusters in 
the Edge sample with PSPC images and published (Tx) 
with uncertainties. Roughly half of these temperature s 
are from Einstein MPC observations ( David et al. 19931) , 
but where possible we substitute more accurate Ginqa. 



Day ei 



ASCA, o r ROSAT PSPC temperatures (GmQa- 
al. 19|9ll lAUcn ct al. 1992|, [Johnstone et al. 1992|, IDavic 
et al. 
1997 
1996 



" , [Hughes et a . 1993- Arnaud fc Evrarc 

ASCA~ \Kenriksen fc Markcvitcli 1996[ [Markevitch 



Matsuzawa ct al. 1996| , Tamura ct al. 19961, iMarkc- 



vitch Vikhlinin 1997| ; PSPC— pavid, Jones fc For ~ 
1996^ We exclude the Virgo cluster from our analysis be- 
cause bright, cluster emission extends beyond the 2° PSPC 
field of view. Thus, our final sample contains 41 clusters. 

The observational data for this sample are displayed in 
Table [|; the cluster name, effective PSPC exposure time 
texp, background brightness Ib, emission weighted mean 
temperature (Tx) , and designation as cooling flow cluster 
are listed. We designate those clu sters with central cooling 
times significantly below 10 Gyr ( Edge, Stewart fc Fabian 
19920as cooling flow clusters. 

Figure Q contains a plot of the 41 clusters in log Ri and 
log (Tx) for I = 1.93 x 10"^ cts/s/arcmin^. A hnear fit to 
the sample (minimizing the sum of the orthogonal residu- 
als from the relation) yields 



logi?j = (0.93 ±0.11) log 



(Tx) 
6keV 



(0.08 ±0.01). (2-2) 



The rms scatter in log Ri about this relation is 
SlogRraw = 0.084. The coefficient uncertainties reflect 
the half width of the 68% confidence region determined by 
bootstrap resampling. 



We analyze the ST relation in the PSPC sample over 
a factor of 8 in surface brightness (see Table || and §3.2); 
varying / primarily changes the zeropoint b of the ST re- 
lation. The integrity of the relation is preserved until / 
reaches values comparable to the peak surface brightnesses 
of some clusters. 

The uncertainties in {Tx) are large enough to con- 
tribute to the observed scatter. We estimate the intrinsic 
scatter in S log Rj by subtracting the temperature con- 
tribution in quadrature. We find A{S log Rj) = 0.93 * 
RMSiar/Tx/ In 10) = 0.054. Thus, the intrinsic scatter 
is approximately SlogRmt = 0.064. Fig. || gives the im- 
pression that the scatter increases at lower temperature; 
fitting only those 30 clusters with {Tx) > 4 keV leads to 
an ST relation logi?/ = (0.82 ± 0.14) \og{{Tx) /6keV) - 
(0.07 ± 0.02), with raw/intrinsic scatter of 0.067/0.041. 
Thus, the slope of the ST relation for hot clusters is slightly 
shallower than that for the whole population, and the scat- 
ter is smaller by 50%. 

3. SIMULATIONS 

We further examine the ST relation using a set of 48 
high resolution N-body plus gas d ynamics simulations, t o 
be described in detail elsewhere (Mohr & Evrard 1997). 
In summary, the cluster simulations are carried out within 
four different cold dark matter (CDM) dominated cos- 
mologies (1) SCDM (f7 = 1, (78 = 0.6, h = 0.5), (2) OCDM 
(f7o = 0.3, 0-8 = 1.0, h ^ 0.8), (3) LCDM (f^o = 0.3, 
Ao = 0.7, 0-8 = 1-0, h = 0.8), and (4) ZCDM {9. = 1, 
as — 1-0, h — 0.5) with power spect rum consistent wit h 
a CDM r = 0.24 transfer function ( pavis et al. 1985[ ). 
Here Hq = lOO/i km/s/Mpc and cts is the power spectrum 
normalization in 8/i^^ Mpc spheres. The baryon density 
is a fixed fraction of the total ilf, = 0.2ilo- Within each of 
these models, we use two 128^^ N-body only simulations 
of cubic regions with scale ~ 400 Mpc to determine sites 
of cluster formation. Within these initial runs the virial 
regions of clusters with Coma-like masses of 10^^ M0 con- 
tain ~10'^ particles. 

Using the N-body results for each cosmological model, 
we choose clusters for additional study. We zoom in on 
these clusters, resimulating them at higher resolution with 
gas dynamics and gravity on a 64'^ N-body grid. The 
simulation scheme is P3MSPH (Evrard 1988), and radia- 
tive cooling is ignored. These clusters have masses which 
varying by an order of magnitude. The scale of the sim- 
ulated region surrounding each cluster are in the range 

1/3 

50-100 Mpc, and vary as M^^;^, where Mhaio is approx- 
imately the mass enclosed within the present epoch turn 
around radius. Thus, the 48 simulated clusters in our final 
sample have similar fractional mass resolution; the spatial 
resolution varies from 125-250 kpc. We create X-ray im- 
ages and temperature maps for further analysis following 
procedures described in Evrard (1990). 

Fig. ||A is a plot of the ST relation for the simulated 
clusters (with / = 1.62 x lO'^^/ifo cts/s/arcmin^). Clus- 
ters from the four cosmological models are plotted with 
different point styles, each appears three times from or- 
thogonal projections. The best fit relation for the SCDM 
model is logi?/ = 0.70 log ((Tx) /6keV) + 0.119 and the 
scatter around this fit is 5logR = 0.042; Table | contains 
the relations and scatter for all models; the level of scat- 
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ter is typically somewhat smaller than the observed value. 
Only ZCDM clusters have a higher scatter than our esti- 
mate of the intrinsic scatter in the observed sample; the 
higher scatter is caused by a few ongoing major mergers 
at high {Tx). 

The scatter around the ST relation in the simulated pop- 
ulations is due to projections of clusters along the line 
of sight, small equilibrium departures, and small (~ 5%) 
intrinsic fractional ficM variations which exist because 
of different formation histories. On average, these clus- 
ters reached half of their present epoch mass at look-back 
epochs O.aato: 0.35to: 0A2to: OAAto: for Z:S:0:LCMD, 
where to is the age of the universe in each model. The 
S, O and LCDM clusters are morphologically similar to 
observed clusters, and the ZCDM clusters exhibit some- 
what more evidence for recent merg ers than the other thre e 
models and the observed clusters (Mohr fc Evrard 1997). 



4. ANALYSIS 

The existence of the ST relation in this large, approxi- 
mately flux limited cluster sample demonstrates that the 
ICM structure is a well behaved function of (Tx). Af- 
ter correcting for measurement uncertainties, the scatter 
of the clusters about the ST relation is similar to the 
scatte r of bright ellipticals around the fundamental plane 
' e.g., J0rgensen, Franx fc Kjg;rgaard 1996| , Mohr & Weg 



ner 1997). To further understand this relation, we examine 
sources of scatter about it, and then we use the scatter to 
place upper limits on variations in the ICM mass fraction 
fiCM ■ We also examine the ST relation slope, which pro- 
vides information about the physical processes affecting 
the ICM structure. 

4.1. Scatter About the ST Relation 

There are several sources of scatter in the ST relation be- 
sides the measurement uncertainties discussed in §2. These 
include: (1) variations in galactic Nh absorption, (2) de- 
partures from equilibrium, (3) projections of physically 
unassociated clusters along the line of sight, and (4) dark 
matter and ICM structural variations at constant (Tx). 

The contribution from variations in galactic Nh is in- 
significant. We use PROS tasks to calculate surface bright- 
ness correction factors for galactic absorption at 1 keV for 
our cluster sample. We find that these correction factors 
have a mean of 1.17 and an rms variation of 6.8% about 
this mean value for the 41 clusters. Typical surfac e bright 



ness profiles are rather steep I(R ) on R ^ (e.g. lones fc 



Forman 1984, Mohr et al. 1995), so 5\ogR = Mlog/ 



and the variation in galactic absorption along the hues of 
sight to these clusters is expected to contribute scatter at 
the level of ~2.3%. Galactic absorption is approximately 
45% more effective at 0.5 keV, the low energy end of our 
bandwidth (and, of course, less effective at 2.0 keV), so 
the maximal contribution to the ST relation scatter from 
variations in galactic Nh is ~3.3%. 

Departures from equilibrium contribute significantly to 
the scatter. Two clusters of similar mass about to merge 
in the plane of the sky will be 40% too large in Rj for their 
(Tx). The scatter about the ST relations exhibited by the 
simulations described in §3 is an indicator of the expected 
contribution from equilibrium departures and projection 
effects in cluster populations still merging at the present 
epoch. 



Equilibrium dark matter and ICM structural variations 
at a particular temperature will also contribute to the 
ST relation scatter. N-body simulations indicate that 
the dark matter structure of clusters is a regular func- 
tion of mass, and that mild variations among clusters with 
similar mass are attributable to different formation times 
( Navarro, Frenk fc White 1996 ). Again, major mergers 



present exceptions, but do not dominate the local pop- 
ulation. The ICM distribution of clusters with similar 
masses has not been as extensively studied as the dark 
matter, but the experiments shown here and below sug- 
gest that the ICM structure is a well behaved function 
of cluster mass or temperature, even in the presence of 
galactic winds (Metzler & Evrard 1997). Other sources of 
variation in the ICM distribution or gas mass fraction are 
possible. One source is radiative cooling, which leads to 
strong emission excesses or cooling fl ows in cluster centers 
dJones fc Forman 1984| , [Fabian 199| ). 

We examine the effects of cooling flows by dividing our 
cluster sample into those with cooling times {tcooi signif- 
icantly less than 10 Gyr) short enough that cooling flows 
should exist in the core, and those with longer cooling 
times ( [Edge, Stewart fc Fabian 19921 ). We find that the 



best fit ST relation for the 18 cooling flow clusters (see 
Table [^) has a slope m = 0.93 ± 0.16 and a zeropoint 
b = —0.07 ± 0.02; the relation for the 23 non-cooling 
flow clusters has m = 1.00 ± 0.21 and b = -0.09 ± 0.02. 
The raw/intrinsic scatter about these two relationships is 
SlogRi = 0.086/0.071 for the cooling flow clusters and 
S log Ri = 0.083/0.061 for the non-cooling flow clusters. 
These two relations are indistinguishable and the scatter 
about each is comparable; thus, for the present sample of 
41 clusters there is no indication that the ST relation is 
sensitive to the presence of cooling flows. 

4.2. Upper Limit on ICM Mass Fraction Variations 

By assuming the observed scatter in the ST relation is 
caused solely by varying ICM mass fraction ficM at con- 
stant (Tx) , we estimate an upper limit on these variations. 
The surface brightness profiles of both real and simulated 
clusters are generally well fit with the standard, isothermal 
/3-model 



I{R) = AKplR, [1 + (i?/i?c)'] 



-3/3-1-1/2 



(4-1) 



where po is the central ICM density, A is the specific emis- 
sivity, Rc is a core radius, /? is the variabl e slope, and 
A is a constant with some (3 dependence ( Cavaliere fc 
Fusco-Femiano 1978 ). We use this model as a guide to 
understanding the ST relation (but the relationship below 
can be derived independent of the (3 model). Within this 
framework, the ICM mass fraction is 

47rpoi?g/„^-/^-dAA^(l + A^)-3/^/^ 
^^^^^ ^ ^'"'^ 

where Rmr is the virial radius encompassing a fixed density 
contrast. Outside the core, the surface brightness scales 
as 

I{R) ^ C f]cM R-''^' (4-3) 

where the constant C incorporates the ICM structural in- 
formation of equation (4-2). If clusters of a given tempera- 
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ture are structurally similar, then the scatter in the isopho- 
tal radii at constant (Tx) places an upper limit on varia- 
tions in the gas fraction (5 log //cm < [(6/3 — l)/2]d\ogRi. 
This inequality can be broken if structural variations (vari- 
ations in C) are anti-correlated with variations in ficM', 
we assume here that this is not the case. Qualitatively, 
such an anti-correlation would require gas-poor clusters 
to be more centrally concentrated (and vice- versa), a sit- 
uation at odds with traditional expecta tions. 
F or the typical val ue /? = 2/3 (^c 



Formaij 

19841, |Mohr et al. 1995|), 5 log //cm < 1.55 log i?/, and the 



intrinsic scatter in logi?/ limits fractional gas variations 
to Sf/ficM < 22%. The limit drops to Sf/fwAi < 14% 
for the 30 clusters hotter than 4 keV. 

The scatter about the ST relation in the simulated 
clusters includes a contribution from variations in //cm 
within the virial radius. This contribution is rather small: 
df / ficM — 4% within rsoo, the radius within which the 
mean cluster density is 500 times the critical density. 
Thus, other effects such as departures from equilibrium 
and projections dominate the scatter. Assuming these ef- 
fects exist to a similar degree in the observations, then the 
similar scatter in the simulations and the observed clusters 
with (Tx) > 4 keV indicates that gas fraction variations 
within hot clusters are very small, perhaps as small as 4%. 
This argument does not apply for cooler clusters, which 



exhib: t larger scatter about the ST relation. 



4.3. Slope of the ST Relation 

The steeper observed slope m = 0.93 ± 0.11, compared 
to 0.61 < m < 0.81 in the simulations, suggests structural 
differences in real clusters compared to this set of mod- 
els. Radiative cooling is not included in the simulations, 
but the similarity of the ST relation in clusters with and 
without cooling flows indicates this is not the cause of the 
slope differences. Possible explanations include systematic 
trends in the overall gas fraction ficM and/or systematic 
variations in the ICM distribution (as measured by /3) with 
(Tx). 

The observed cluster sample provides evidence that shal- 
lower surface brightness profiles lead to steeper ST rela- 
tions. Table || lists the best fit slope m and zeropoint b 
for the observed clusters over a factor of 8 variation in /. 
Along with the expected zeropoint decrease with brighten- 
ing /, there is a tendency for a steeper slope. This steep- 
ening is related to the shallower surface brightness profiles 
(smaller effective /3) at brighter isophotes. Simulations 
from the four cosmological models provide additional ev- 
idence that shallower surface brightness profiles generally 
lead to steeper ST relation slopes; the mean f3's of the 
Z:S:L:OCDM are are 0.67:0.73:0.77:0.81 and the m's are 
0.81:0.70:0.72:0.61. Thus, the steeper observed ST rela- 
tion may simply be an indication of a more extended ICM 
distribution in the observations than in the simulations. 

Galactic winds are one physical process expected to in- 
troduce systematic ICM structural variations. In such a 
model, star formation within galaxies results in supernovae 
(SNe) driven winds which expell gas, and the expelled gas 
has the orbital kinetic energy of its parent galaxy plus 
the energy imparted by the SNe. This non-gravitational 
source of energy is expected to be more important in low 
(Tx) clusters where the orbital energies are smaller. We 
examine the effects of galactic winds on the ST relation 



using an ensemble of clusters from simulations which in- 
clude galaxy feedback (Metzler & Evrard 1994; 1997). The 
simulated cluster ensemble contains two populations; one 
population of 18 clusters evolved consistent with an SCDM 
model and another population with the same initial condi- 
tions and cosmological parameters, but with galaxy feed- 
back modeled by discrete gas ejection totalling half the 
initial galaxy mass between z = 4.5 and the present. This 
extreme ejection model is intended to estimate the maxi- 
mal effect of galaxy feedback. 

Fig. ||B is a plot of the effects of galaxy feedback on the 
ST relation. The slope of the best fit ST relation for the 
clusters with no feedback (dotted line) is to = 0.69, con- 
sistent with the SCDM slope in our simulations (see Table 
|2|). The slope for clusters simulated with ejection (sohd 
line) is TO = 0.99, somewhat steeper than the slope of our 
PSPC sample. This demonstrates that galaxy feedback 
can cause the structural changes required to match the 
slope of the observed ST relation. Structural changes in- 
clude shallower surface brightness profiles (smaller 0) and 
lower overall gas fractions with decreasing {Tx), as dis- 
cussed in Metzler & Evrard (1997). Fig. ||B makes clear 
that, as expected, ejection has the greatest effect on the 
lowest mass clusters. 

The tendency for low (Tx) clusters to have shallower 



surfa ce brightness profiles is well known (e.g., Mohr et al 
I995J); our sample of 41 clusters exhibits this general trend. 
Fig. g contains a plot of the slope of the surface brightness 
profile Pef f measured over the region of the cluster used 
to cal culate Rj . The S pearman rank correlation coeffi- 
cient ( [Press et al. 1992| ) for (Tx) and /3eff is = 0.54 

Peff is a two 
is measured 



with a correlation significance of 99.97%. 
dimensional generalization of /3, and like Rj 
without azimuthal averaging; moreover, Peff reflects the 
local slope of the surface brightness profile and so is in- 
dependent of the cluster core radius. Briefly, to calculate 
the values in Fig. |^ we determine the area Ai enclosed by 
the isophote /; = 1.29 x 10^'^ cts/s/arcmin^ (chosen to be 
somewhat fainter than the / used to calculate Rj ), and 
then calculate the isophote Ih which encloses the area Ah 
which is chosen to be Ai/2. This provides two isophotes 
and two areas, and so we then calculate the effective slope 
of the surface brightness profile between /; and Ih 



Peff 



1 / log(/^//0 
3 \logiAi/Ah) 



(4-4) 



This form follows directly from the expression I{R) oc 
^-6/3,^^+1 rpj^g ^^g^^ falling P^ff with (Tx) is qualita- 
tively consistent with the effects of galaxy feedback. 

5. DISCUSSION 

We use observations of the 41 members of an X-ray flux 
limited cluster sample with PSPC observations and pub- 
hshed (Tx) (except Virgo) to demonstrate that nearby 
clusters conform to a tight relation between X-ray isopho- 
tal size and cluster temperature. The intrinsic scatter in 
size at fixed temperature is only 15% (10% for clusters 
with (Tx) > 4 keV). The existence of the ST relation in- 
dicates that the ICM structure outside the core regions is 
a well behaved function of (Tx) and suggests a tight cor- 
relation between (Tx) and M^r as seen in our simulations 
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and others (Evrard, Metzler fc Navarro 1996, Schindlei 



1996 1, iRocttigcr, Burns fc Loken 199^ ) 

The ST relation is significantly tighter than the Lx — 
{Tx) relation for this cluster sample. The scatter (cor- 
rected for t emperature uncer tainties) in 2-10 keV X-ray 
luminosity (David et al. 1993) around the best fit relation 
for these 41 clusters is 5Lx / Lx ~ 52%. (We use the same 
set of temperature measurements in both cases.) The fac- 
tor of 3.5 larger scatter in the Lx — {Tx) relati on refiects 
the sensitivity of Lx to cluster core properties (Fabian et 
al. 19^). When restricted to weak cooling flow clusters — 
those with inferred cooling flow rates below 100 yr~^ 
— the Lx — {Tx) relation displays a much smaller disper- 
sion, and inferred gas fraction variations from a sample of 
24 clusters are consistent with the limits set by this analy- 
sis (Arnaud & Evrard 1997). The narrow scatter in the ST 
relation indicates that structural regularity exists outside 
the core at observationally viable surface brightness levels. 

Our cluster sample provides no evidence that cooling 
flows significantly affect the ST relation. The 13 clus- 
ters with central cooling times significantly below 10 Gyr 
exhibit a best fit slope, zeropoint, and scatter which are 
statistically indistinguishable from those for the 18 other 
cluste rs. The zeropoint differences between these two ST 



relations indicate that {Tx) in cooling flow clusters is bi- 
ased low by ^ 5 ± 7%, consisten t with previous estimates 
of -10-20% ( [Fabian et al. 1994] ). 

The scatter around the ST relation is comparable to 
that of the half light radius of elliptical galaxies around 
the f undarncntal plane (e.g. J0rgensen, Franx fc Kj33r-| 
gaard 199£, Mohr fc Wegner 19971 ); an obvious use of the 
ST relation is as a distance indicator. Knowledge of {Tx) 
and an apparent isophotal size predicts the cluster distance 
with an uncertainty limited by the scatter in Ri 10% 
for hot clusters). We are currently evaluating the promise 
of this relation as a distance indicator at intermediate red- 
shift. 

The scatter about the best fit ST relation for observed 



clusters provides an upper limit on the rms variation 
in cluster ICM mass fraction at a common temperature 
{Tx)- Ignoring possible correlations between cluster struc- 
tural changes and ficM variations, we find upper limits 
Sf/ficM < 22% (< 14% for {Tx) > 4 keV). 

We use 48 cluster simulations within 4 cosmological 
models to demonstrate that cluster populations experi- 
encing growth at the present epoch exhibit ST relations 
with scatter similar to that observed above {Tx) > 4 keV: 
SRj/Ri ~ 10%. Only a small contribution to the scatter 
in the simulations is from ficM variations, suggesting that 
the ICM mass fraction in hot clusters may be limited to 
< 5% variation within their virial regions. Implications 
for hierarchical models of galaxy formation remain to be 
explored. 

The slope of the ST relation for the PSPC clusters is 
steeper than the slope in our simulations. We use numer- 
ical simulations (Metzler fc Evrard 1994, 1997) to show 
that galaxy feedback introduces the kinds of ICM struc- 
tural changes required to steepen the ST relation. The 
changes include shallower surface brightness profiles and 
lower ficM with decreasing {Tx)- A decrease in /3 with 
(Tx) ( see Fig. | 3[) and a reduction of f i cM in low {Tx) sys- 
tems (David, Jo _ nes fc Forman 1995 , dell' Antonio, Gellei 



fc Fabricant 1995) have been observed. Galactic feedback 



is a promising mechanism for creating the observed, struc- 
tural differences between low and high (Tx) clusters. 
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Fig. 1. — The X-ray ST relation for 41 rliist ers. The cluster redshift is used to calculate Rj , and / = 1.93 X 10~^ cts/s/arcmin^; error 
bars are Icr. The best fit ST relation (see Eqn 2-2) is plotted, and the rms scatter around this relation is (51ogi?/ = 0.084. The intrinsic 
scatter about the relation (removing contribution trom {Tx) uncertainties) is SlogRj = 0.064 (15%); for 30 clusters with (Tx) > 4 keV, the 
scatter is <51ogR/ = 0.041 (10%). 
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Table 2 

Fit Results for Size-Temperature Relation 
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Fig. 2. — The ST relations for simulated galaxy clusters. Our 48 simulations appear on the left (A), with the points coded according 
to cosmological model. The best fit relation for the SCDM model is plotted (see Table H), and the rms scatter around this relation is 
& log K = 0.042. The scatter for these simulated clusters (which have formed recently) is similar to that for observed clusters, but the slope of 
the relation is shallower. On the right (B) we demonstrate the effects of galaxy feedback on the ST relation with simulations from Metzler & 
Evrard (1994, 1997). The best fit slope for the clusters simulated with feedback (m = 0.99) is much steeper than the slope of those simulated 
without feedback (m = 0.69). 
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Fig. 3. — Slopes of the X-ray surface brightness profiles are plotted versus cluster temperature {Tx)- The trend for cooler clusters to have 
shallower profiles is qualitatively consistent with the effects of galaxy feedback. j is calculated independent of the cluster core radius and 
requires no aaimuthal averaging (see §4.3 for details). 



